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(54) Abstract Title 

Gradually deforming a stochastic model of a heterogeneous medium eg underground zone 

(57) Method for gradually deforming, totally or locally, a Gaussian or similar (i.e. Gaussian-based) stochastic 
model of a heterogeneous medium such as an underground zone, constrained by a series of parameters 
relative to the structure of the medium, comprises drawing a number p, at least equal to two, of independent 
realizations of at least part of the model, and linear combination of these p realizations with p coefficients such 
that the sum of the squares of the coefficients is equal to 1 . This linear combination constitutes a new 
realization of the stochastic model and it gradually deforms when the p coefficients are gradually modified. 
TTie method can comprise several iterative gradual deformation stages, with combination at each stage of a 
composite realization obtained at the previous stage with q new independent realizations. The method makes 
it possible to gradually deform realizations of a model representative of the medium while modifying the 
statistical parameters relative to the structure of the medium. The method also allows gradual individual 
deformations of various parts of the model while preserving continuity between these parts. The method can 
be applied for construction of stochastic reservoir models constrained by non-linear data (data linked with the 
flow of fluids). 
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At least one drawing originally filed was informal and the print reproduced here is taken from a later filed formal copy. 



GB0gg342736[httM^ 



1/6 




Fig.l 



GB002342736fhttp://www.q fltthppatent.t«m/LoQin.doa /$cbmnd 




GB0Q234273 6[htto /A«u^natthe »»atent.com/Login4O9 g.^M]didgft/ 



rage h oi 



3/6 




le^Ol le^02 IrKii le*04 le+05 le*06 



Fig.3b 



r.Rm^-^d97aRthWpr/A««w.»etth6 Patentcom/Loq in.dogg^n^ 



4/6 




le-01 H 1 1 1 1 

le+01 le'»'02 le-f04 le^-OS le+06 



Fig.4b 



GB002342736 r http:/Awww.qetthepatent.com/LoQin.doq/$cbmndidqe/Fetch/GB0023427 part=maintooiDar=Dottom i 



Pa g e fa oT 20 



516 




Xe-01 I I I I I I I I I I j I I 1 I { I I I 



-1.0 -0.5 0.0 0.5 1.0 



Fig.5b 



GB002342736fh»p:/AAMw.o etlheoatentcom/Lo Qin.doa/$cbmndidq e/Fetch/GB002342736.cpcm^ 



Pa g e ( or /s^ 




GB00234273P rhttp://v^ww.qetthepatent.com/Login.dog/$cbmndidge/Fetch/G^ Pa g e 8 or 

2342736 

METHOD FOR GRADUALLY DEFORMING A STOCHASTIC MODEL OF A 
HETEROGENEOUS MEDIUM SUCH AS AN UNDERGROUND ZONE 
The present invention relates to a method for gradually deforming, partly or totally, 
realizations or representations of a heterogeneous medium such as an underground zone, 
which all are the expression of a Gaussian or similar type stochastic model, i.e. 
5 including an underlying Gaussian model. 

The method according to the invention finds applications notably in the construction 
of a stochastic model of an underground formation, constrained by non-linear data. 

BACKGROUND OF THE INVENTION 

Examples of use of Gaussian models for modelling the subsoil structure are 
10 described by : 

- Joumel, A.G. and Huijbregts Ch J. : "Mining Geostatics", Academic Press, 1978, or 

- Matheron, G. et al. : "Conditional Simulation of the Geometry of Fluviodeltaic 
Reservoirs", paper SPE 16753; SPE Annual Technical Conference and Exhibition, Las 
Vegas, 1987. 

15 Constraint of stochastic models by a series of non-linear data can be considered to 
be an optimization problem with definition of an objective function assessing the 
accordance between data measured in the medium to be modelled and corresponcUng 
responses of the stochastic model, and minimization of this function. A good 
deformation method must respect the following three characteristics : a) the deformed 

20 model must be a realization of the stochastic model; b) it must produce a regular 
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objective function that can be processed by an efficient optimization algorithm; and c) 
the deformation must be such that the solution space can be entirely covered. 

The method according to which values are exchanged at two points of a realization 
of the stochastic model, i.e. a realization or representation of the medium studied that is 
5 the expression of the model, can be cited as an example of known deformation method. 
This method is used in the technique referred to as "simulated annealing" known to 
specialists. An example thereof is described by : 

- Deutch, C. : "Conditioning Reservoir Models to Well Test Information", in Soares, 
A. (ed.), Geostatistics Troia*92, pp.505-518, Kluwer Academic Pub., 1993. 

10 As an arbitrary exchange breaks the spatial continuity of a stochastic model, it is 
necessary to include a variogram in the objective function, which makes optimization 
very tedious. 

Another known method is the method referred to as pilot point method, which is 
described for example by : 
15 - de Marsily. G. et al. : "Interpretation of Interference Tests in a Well Held using 
Geostatistical Techniques to Fit the Permeability Distribution in a Reservoir Model" in 
Verly. G. et al. (ed.), Geostatistics for Natural Resources Characterization, Part 2, 831- 
849, D. Reidel Pub. Co, 1984. 

This method essentially consists in selecting a certain number of points of an initial 
20 realization (representation), in calculating the derivatives (sensitivity coefficients) of the 
objective function in relation to the values at these points, and then in modifying the 
values at pilot points in order to take account of these sensitivity coefficients. A new 
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realization is formed by the known method referred to as conditional kriging. It may 
follow therefrom that a change in the values at the pilot points, by using the sensitivity 
coefficients, leads to an undue change in the variogram even if conditional kriging is 
used afterwards, in particular when the number of pilot points becomes great. 

5 SUMMARY OF THE INVENTION 

The method accoiding to the invention allows to carry out at least partial gradual 
deformation of the realization (or representation) of a Gaussian or similar type 
stochastic model of a heterogeneous medium such as an underground zone, constrained 
by a series of parameters relative to the structure of the medium. It is characterized in 
10 that it comprises drawing a number p at least equal to two of independent realizations of 
at least part of the selected medium model from all the possible realizations, and linear 
combination of these p realizations with p coefficients such that the sum of tfieir squares 
is equal to 1. This linear combination constitutes a new realization of the stochastic 
model and it gradually deforms when the p coefficients are gradually modified. 

15 More generally, the method can comprise several iterative stages of gradual 
deformation, with combination at each stage of a composite realization obtained at the 
previous stage with q new independent realizations drawn from all the realizations. 

The method allows gradual deformation of a model representative of the medium 
while modifying the statistical parameters relative to the medium structure. 

20 The method also allows to perform individual gradual deformations of various parts 
of the model while preserving continuity between these pans (i.e., for example, 
reproduction of a variogram linked with this model). 
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The coefficients of the linear combination preferably are trigonometric functions. 

A Gaussian model or a model similar to the Gaussian type can be used as the 
model : lognormal model, truncated Gaussian model, etc. 

In relation to the prior art, the method according to the invention is founded on a 
5 more solid theoretical basis for stochastic model deformation. It also makes it possible 
both to modify the parameters (often a priori uncertain values) of variograms and to 
deform a realization of a stochastic model. Fmally, thanks to the possibility of 
individual deformation of various parts of the model, the method according to the 
invention gives greater flexibility and efficiency to the operator for adjustment of the 
10 stochastic model to the field reality. The method according to die invention thus allows 
to esublish a connection between stochastic model adjustment and deterministic model 
adjustment by zoning conventionally used by reservoir engineers. 

BRIEF DESCRIPTION OF THE FIGURES 

Other features and advantages of the method according to the invention will be clear 
15 from reading the description hereafter of a non limitative example, witii reference to the 
accompanying drawings wherein : 

. Figure 1 shows an example of an iterative optimization pattern according to die 
method, 

. Figure 2a shows an example of lognormal permeability field tiiat constitutes a 
20 reference model, 

. Figure 2b shows the pressure and logarithmic derivative curves obtained from a 
numerical simulation of a well test in the reference model. 



I 
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- Figure 3a shows an example of a lognormal permeability field associated with a first 
realization of the model, 

- Figure 3b shows a comparison between the pressure curves associated with 
realization 1 and the reference model respectively, 

5 - Figure 4a shows an example of a lognormal permeability field associated with a 
second realization of the model, 

- Figure 4b shows a comparison between the pressure curves associated with 
realization 2 and the reference model respectively, 

• Figures 5a and 5b respectively show the variation curves of the objective function as 
10 a function of parameter p in a linear scale and in a logarithmic scale, 

- Figure 6a shows an example of a lognormal permeability field associated with a 
realization constrained by application of the deformation method, 

- Figure 6b shows a comparison between the pressure curves of the constrained 
realization and of the reference model. 

15 DETAILED DESCRIPTION OF THE INVENTION 

General points 

The method according to the invention allows to gradually deform a stochastic 
model while respecting its statistical characteristics. This method is operational within a 
Gaussian framework. Combined with an optimization method, this method constitutes 
20 an efficient tool for calibrating an image of the medium to be modelled to non-linear 
data. 
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Gradual deformation at the set statistical'parameters 
Consider a series of possible realizations Rj of a stochastic model of the medium to 
be modelled defined by a certain set of representative statistical parameters (average, 
variance, variogram, etc.). These various realizations are obtained from drawing a 
5 random germ; two successive values of this germ can lead to very different realizations. 

Gradual deformation of the model is obtained according to the method by a series of 
iterative combinations of realizations R, of the model. As shown in Figure 1, the basic 
combination consists, from two realizations R„ Ri obtained by drawing, in forming a 
realization Rci meeting the relation as follows : 

10 Rci = a,.Ri + ai.Ri (1) 

with the condition : 
a,-+ai' = I (2) 

This combined realization gradually deforms when coefficients a, and Oa are 
gradually modified. 

On the next iteration, another reaUzation Rj is drawn and combined according to the 
15 same rule with tiie previous combined realization R«|. 

Rji = aj.R«, + cx».R3 

f • 

The new coefficients of the combination also meet die condition : 



fageiJOLsy 
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If the realization Rj drawn leads to a combined realization that does not sufficiently 
match the measurements, one or more successive draws are performed until a 
satisfactory realization is obtained. 

The previous combination pattern can be readily extended to any number of 

5 realizations. At the first combination stage, p different realizations are drawn and 

combination is performed according to the relation as follows : 

R^Y,(X,R, (3) 

with the condition : 

I«f = l (4) 

i 

At each of the next stages, the previous combined realization and q other 
realizations obtained by drawing are combined and combination thereof is achieved 
10 according to the same rule. 

According to an advantageous implementation mode, two realizations Ri, R2 

connected by the relation : 

/? = /?, cos(p;r) + jRj sin(p;r) (5) 

that meets condition (2) are combined at each stage. ^ 

Minimization of an objective function allows to determine coefficient p. It is a 
15 mono-parameter optimization problem. 

In the more general case where p+1 realizations are used, we have a problena of 
optimization of the p parameters Pi 0=1,2.... ,p). and coefficients Oj (i=0,1.2„..,p) are 
written in the following form : 
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a, = sin( A;r) flcosiPjH) (i = 1. p - D 



a^=sin(p/r) 



5 that meet coiidition (4). 

Gradual deformation by modifying tlie statistical parameters 

If no precise values of the parameters of a stochastic model are available, these 
parameters must also be identified by inversion. The method according to the invention 
allows gradual deformation of a reaUzation of a stochastic model by modifying the 
10 parameters of the variogram. To this end. a realization is written in the following form : 

R = L[Y1 

where Y is a white noise (realization of a Gaussian field without correlation), and L an 
operator determined by the variogram and allowing to turn Y into a realization R 
respecting this variogram. Tlie techniques to be implemented here are well known to 
15 specialists. One of them is described by : 

. Oliver. D.S. : "Moving Averages for Gaussian Simulation in Two or Tthree 
Dimensions", Math. Geology. Vol. 27. No.8. 1995. 

We can thus deform Y as before while modifying the variogram operator L : 

R = L[Y1 = L[Y, cos(pn) + Yi sin(p7t)l 



I 



GB0 023427^ [http/AAAAAA/gBtthftpatfint rnm/LQqin.doa /$cbmndidqe/Fetch/GB00234^ Page16of20 

9 

where Y| and Yi are two white noises independent of each other. 
Gradual deformation by parts 

In order to improve efficiency of the calibration of a model to non-linear data 
distributed in the medium to be modelled, one may be led to divide the medium into 

5 several parts and to perform calibration part by part. The method according to the 
invention allows to cany out individual gradual deformations of various parts of the 
medium to be modelled while preserving continuity between these parts. Suppose that 
the medium is divided in n parts with which n white noises Yi, Y:,..., Ya are associated. 
These white noises can be deformed individually by means of the previous algorithm 

10 and correlated by correlation operator L. 

K cos(p,;r) + sinCpjTr) 
};cos(p,7r) + y;sin(p^;r) 

15 where Yi*and Y " (for i=l,2,. , .,n) are two white noises independent of each odier. 
Extension to the other models 

The gradual deformation method described above can be used in the case of models 
transformed from a Gaussian model such as the lognormal model and the truncated 
Gaussian model. 



/? = L[y]=L 
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Validation example 

The method according to die invention is illustrated here by a synthetic study of the 
constraint of a permeability field by pressure dau from a well test. A synthetic 
permeability field (Fig.2a) constmcted according to a lognormal stochastic model is 
5 used as a reference model. A numerical well test is carried out on this field and the 
result is a pressure curve and its logarithmic derivative as shown in Fig.2b. What has to 
be done is to construct a realization constrained by this pressure curve odier than the 
reference image. 

Figures 3a and 4a show two independent realizations of the same model. It can be 
10 observed that the pressure curves of the well tests on these two realizations are very far 
from those of the reference model (Figs.3b and 4b). A chain of realizations has been 
constructed from these two realizations according to the rule given by formula (4). 
Figures 5a and 5b show the variation of the objective function as a function of 
parameter p (that ranges between -1 and 1). The two initial realizations respectively 
15 correspond to p = Oandp = 0.5. Itcan be noticedthatdte objective fiinction reaches its 

minimum value with p = -0.3 or thereabouts. 3 to 4 iterations are sufficient to minimize 
this function when using the gradient method. Figure 6 shows the realization 
corresponding to p = -0.275 obtained by optimization and it can be seen that calibration 
of this realization to the well test is very satisfactory. 
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CLAIMS 

1) A method for gradually deforming at least partly realizations of a Gaussian or 
similar type stochastic model of a heterogeneous medium such as an underground zone, 
constrained by a series of parameters relative to the medium structure, characterized in 
5 that it comprises drawing a number p, at least equal to two, of independent realizations 
(Ri) of at least part of the selected model of the medium from all die possible 
realizations, and at least one stage of gradual deformation of the model by linear 
combination of the p initial independent realizations, coefficients (Oi) of this 
combination being such that the sum of their squares is equal to 1. 

10 2) A method as claimed in claim 1, characterized in that it comprises several 
iterative gradual deformation stages, with combination at each stage of a composite 
realization obtained at the previous stage with q new independent realizations drawn 
from all the realizations. 

3) A method as claimed in claim 1 or 2, characterized in that it comprises gradual 
15 deformation of a model representative of the medium while modifying the statistical 

parameters relative to the medium structure. 

4) A method as claimed in claim 1 or 2, characterized in that it comprises 
individual gradual defonnations of different parts of a model representative of the 
medium while preserving continuity between these parts. 

20 5) A method as claimed in any one of the previous claims, characterized in that 
coefficients (cti) of the linear combination are trigonometric functions. 
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6) A method as claimed in ai^ one of the previous claims characterized in that a 
5 model similar to the Gaussian type such as a lognormal modd or a truncated Gaussian 
model is used. 



7) A method for gradually drfbrming at least partly realizations of Gausaan or 
similar type stodiastic modd of a heterogeneous medium such as underground zone 
substantially as hereinbefore described with reference to the acconq)anying drawings. 
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